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ABSTRACT 



Classically time is kept fixed for infinitesimal variations in problems in me- 
chanics. Apparently, there appears to be no mathematical justification in the 
literature for this standard procedure. This can be explained canonically by 
unveiling the intrinsic mathematical structure of time in Lagrangian mechan- 
ics. Moreover, this structure also offers a general method to deal with inertial 
forces. 



Contents 

Introduction S 

0. Notations and definitions H 

1. Structure of a second order differential equation relative to a metric [s] 

2. Newton-Lagrange mechanics of a free system [lo| 

3. Constrained systems [ll 

4. The class of time and the calculus of variations [16 

5. Time constraints |20 

6. Modification caused by a time constraint in a free system |22 

7. Linear constraints depending on time |26 

8. Reference frames and inertial forces |28 
Acknowledgements. |35 
References 35 



2000 Mathematics Subject Classification. 70A05, 70H03, 70H45, 70F20, 70F25, 37J05, 37J60. 
Key words and phrases. Lagrangian mechanics, time. 



2 



J. MUNOZ DIAZ 



Introduction 

This article studies two related questions in Lagrangian mechanics: the mathema- 
tical nature of time and the one of reference frames and inertial forces. 

An isolated mechanical system without constraints corresponds to a Riemannian 
manifold (M, T2) (the configuration space) and a 1-form a on TM (the work form). 
The general abstract version of Newton law postulates the existence of an associated 
tangent field D on TM, that is, a second order differential equation, which is the 
motion law for [M^T^^oi). There is no function f on M that could be used as 
parameter for all solutions of D. Classically a new function t is added to the system. 
This is "the time", which is measured by some "clock" out of the system. Then M 
is replaced by R x M and to the collection of differential equations, t = 1 is added. 
When constraints (holonomic or not) are introduced in the extended manifold R x M, 
the classical procedure to derive the equations of motion is to apply the principle of 
virtual works using infinitesimal displacements in which the time remains unchanged. 
This is usually done without justification, or providing only "physical" arguments; 
see e.g. 0, pag. 3, [H, pag. 48-49, 0, pag. 522, [6], pag. 65, [8], pag. 215. 

The addition of an external time to the configuration space is somewhat artificial. 
In an isolated system there is no "time" function, but there is a canonical "class of 
time" consisting of all the 1-forms on TM, out of the zero section, that contracted 
with second order differential equations give 1. When these forms are specialized to 
each solution of a second order differential equation their primitive functions serve as 
time parameter. 

The admissible infinitesimal variations in mechanics are described as those tangent 
fields 6 on TM which are infinitesimal contact transformations, projecting to M, and 
preserving the class of time. When properly formulated, this last condition appears 
to be equivalent to the classical commutation formula d o S = S o d. 

Hypersurfaces on TM equipped with local time functions are "time constraints" 
considered in Section 5. When free mechanical system with such constraint are consid- 
ered, the equations of motion are modified in a similar way as for ordinary holonomic 
constraints. In a precise sense time constraints can be proved to be deformable to or- 
dinary holonomic constraints. The precise mathematical justification for keeping the 
time fixed in infinitesimal displacements in D'Alembert principle is given by Theorem 

m 

Our presentation of Lagrangian mechanics provides a natural framework for the 
understanding of general inertial forces. The nature of such forces seems to be an 
obscure point in the classical literature. 

Given a manifold isomorphism ip: TZ ^ Ai, where TZ, Ai are provided with respec- 
tive pseu do- Riemannian metrics T2, T2, for each mechanical structure (A1,T2,a) on 
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M, there correspond two mechanical structures {7l,T2,a), {7l,T2,ai) canonically as- 
sociated by ip to the one in A4. Their difference can be read in TZ as the inertial force 
caused in {Ai,T2,ot) by the reference frame </?. This gives a precise mathematical 
meaning to the expression "(/? preserves the equations of motion" that spreads across 
the literature without a proper previous definition. The property of "preserving the 
equations of motion" is interesting only when time constraints are present, as is the 
case for uniparamctric automorphism groups, because for free systems, (p preserves 
the equation of motion only when it is an isometry, and therefore a simple change of 
coordinates. 

Sections 0, 1, 2, 3 have been included to make the article self-contained, although 
the material therein is classical, except for language. 



0. Notations and definitions 

We start with a brief overview of notations and definitions used in the article. 

Let M be a smooth manifold of dimension n, and TM be its tangent bundle. 

Each differential 1-form a on M can be considered as a function on TM, denoted 
by a, which assigns to each ?;„ G T^M the value a(t'a) = (tta- ' a) obtained by duality. 
In particular, a function / G C°°{M) defines the function on TM associated to df 
that wc denote in short by /. This definition also applies to differential forms a on 
TM that are at each point the pull-back of a form on M. In the sequel we call these 
forms horizontal forms. 

The map f ^ f from C°°{M) to C°°{TM) is a derivation of the ring C~(M) taking 
values in the C°°(M)-module C°°{TM). We denote it by d since it is essentially the 
differential. For each horizontal form a, we have a — {a, d) as functions on TM. 

Any derivation 8 of C°°(M) to the C°°(M)-module C°°(TM) can be viewed as a 
field in TM taking values in TM, that is, as a rule assigning to each point Va G T^M 
a tangent vector 5^^ G T^M. More precisely, 5y^f = {6f){va), for each / G C°°(M). 
In particular, the derivation d is the identity vector field in TM, Va Va, since 

dvJ^{df){Va)^{df,Va)^Va{f). 

Using the vector space structure of the fibers of TM we can associate to each 

Va G TaM a tangent vector to T^M at each point as the derivative along Va in TaM. 
Denoting by Va this derivation, we have for / G C°°{M) and a point Wa G TaM: 

Va{f){Wa) = hm = f{Va) = Va{f ) . 

At each Wa G TaM, Va G {TaM) is called the vertical representative of Va G T^M 
and Va the geometric representative of Va. 
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This canonical association between tangent vectors to M at the point a ^ M and 
tangent vectors to the fiber T^M at each one of its points, estabhsh an isomorphism 
between fields on TM valued on TM, and vertical tangent fields on TM. Under this 
isomorphism, the field d corresponds to the the vertical tangent field D on TM such 
that Df = /, for / G C°°(M). This field is the infinitesimal generator for the group 
of homotheties of the fibers of TM. To avoid any confusion with the notation for 
second order differential equations defined below, we will no longer use the capital 
notation D. 

Definition 0.1. (Second Order Differential Equation). A vector field D on 
TM is a second order differential equation when its restriction (as derivation) to the 
subring C°^{M) ofC'^{TM) is d. 

This is equivalent to have ^^{D^^) = Va at each point Va G TaM (where tt: TM 
M denotes the canonical projection). 

Remark 0.2. The difference between two tangent vector fields on TM which are 
second order differential equations, is a vertical vector field. Thus the second order 
differential equations on TM are sections of an affine bundle modeled on the fiber 
bundle over TM of the vertical tangent fields. This last one is isomorphic to the 
bundle of fields on TM taking values on TM. 

Definition 0.3. (Contact System). The contact system onTM is the P faff system 
in TM which consists of all the 1- forms annihilating all the second order differential 
equations. It will he denoted by Q. 

Remark 0.4. The forms in the contact system also annihilate the differences of sec- 
ond order differential equations, i.e. all vertical fields. Therefore, they are horizontal 
forms; each ojy^ G 0.^^ is the pull-back to T*^TM of a form in T*M. Now, a horizontal 
1-form kills a second order differential equation if and only if it kills the field d. Thus 
the contact system on TM consists of the horizontal 1-forms which annihilate d. 



0.5. Local coordinate expressions. 

We take local coordinates {q^, . . . , g") in M and corresponding {q^, . . . , g", g^, . . . , g") 
in TM. We have, using Einstein summation convention, 

^ 9g* 

A vertical field has the expression 
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And the one for a second order differential equation is 

Usually we will denote /* by understanding that it is a given function of the g's 
and g's. 

A local system of generators for the contact system out of the zero section, is 
given by 

q'dq' - q^dq' (z, j = 1, . . . , n) . 



1. Structure of a second order differential equation relative to a 

METRIC 

Let T*M be the cotangent bundle of M and tt : T*M M the canonical projection. 
Recall that the Liouville form 9 on T*M is defined by 9a^ = 7r*{aa) for G T*M. 
Abusing the notation we can write = aa- 

The 2-form uj2 = dO is the natural symplectic form associated to T*M. 

In local coordinates . . . , g") in M, and corresponding (g^, . . . , g",pi, . . . for 
T*M, we have 

9 = Pidq\ UJ2 = dpi A dq^ . 

Let T2 be a (non-degenerate) pseudo-Riemannian metric in M. Then we have an 
isomorphism of vector fiber bundles 

TM T*M 

{ivaT2 is the inner contraction of Va with T2). 

Using the above isomorphism we can transport to TM all structures on T*M. In 
particular, we work with the Liouville form 9 and the symplectic form uj2 transported 
in TM with the same notation. 

From the definitions we have for the Liouville form in TM , at each Va € TaM, 

(1) K = , 

where the form of the right hand side is to be understood pulled-back from M to 
TM. 

Equation ([T]) is the point expression for: 

(2) 9 = t^T2 
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On the other hand, from ([T]) we obtain, for each A G M: 

(3) Oxva = ^^Va 

(equahty as 1-forms in T*M). 

Let V denote the vertical field corresponding to d; V is the infinitesimal generator 
of the group of homotheties in the fibers of TM, and from ([3]) we get 

(4) Lve = e , 

which, by using Cartan's formula (see |2], P- 36) and iyO = {6 being horizontal), 
gives 

(5) iydO = 9 or ivuj2 = . 

Taking the values at each Va G TaM and if we put together ([1]) and ([5]) it results 
the following key lemma. 

Lemma 1.1. If Va ^ TaM and Va is its vertical representative at each point of TaM, 
we have 

(6) iv^uj2 = iv,T2. 

If V is a vertical tangent field on TM and v is the corresponding field on TM taking 
values in TM, then we have 

(7) iv(^2 = ivT2 
( equality of horizontal forms in TM ). 

Definition 1.2. (Kinetic Energy). The function T = ^9 on TM is the kinetic 
energy associated to the metric T2. So, for each Va G TM, we have T{va) = \d{va) = 
^T2{va,Va) or, as a function on TM, T = ^T2{d,d). 



Theorem 1.3. The metric T2 establishes a one-to-one correspondence between second 
order differential equations and horizontal 1-forms in TM. 

The second order differential equation D and the horizontal 1-form a that corre- 
spond to each other are related by 

(8) iDUJ2 + dT + a = . 

Proof. Given a second order differential equation D, we define the 1-form a by ([8]). 
Now we check that a is horizontal. For any vertical field V we must prove that 
{a, V) = 0. Using Cartan's formula. 



{iDu;2, V) = {lode, V) = D{9, V) - V{9, D) - {9, [D, V]) 
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Now we have (6', V) = since 6 is horizontal. Also we have {6, D) = (6, d) = 6, 
and V6 = 2{6,v), where v is the geometric representative of V and using that 6 is 
homogeneous of the second degree in the variables q. Therefore we have 

V{9,D)=2{9,v) . 

Also for / G C°°(M) we have 

[D,V]f=-ViDf) = -Vf = -vf, 
thus the bracket [D, V] is equal to —v up to a vertical field. It follows that 

{6, [D,V]) = -{e,v) . 
And putting all together, we have 

{tDUJ2, V) = 0-2{e,v) + {6, v) = -{6, v) . 

Also 

{dT,V) = VT = V{U) = {e,v) , 

therefore a is horizontal. 

Lemma (11.11) establishes a linear isomorphism, V ^ iv^2i between vertical vector 
fields and horizontal 1-forms. So, given a, adding a suitable vertical vector field V to 
a given second order differential equation Dq we obtain D = Dq + V that corresponds 
to a by dH]). □ 

Definition 1.4. (Geodesic Field). The geodesic field of the metric T2 is the second 
order differential equation, Dq, corresponding to a = 0: 

(9) idg^2 + dT = . 

The projection to M of the curves solution of Dq in TM are the geodesies of T2. 

The geodesic field Dq is chosen as the origin in the affine bundle of second order 
differential equations. With this choice we establish a one-to-one correspondence 
between second order differential equations and vertical tangent fields. 

D< — > D-Dg = V. 

And recalling that to a vertical field V there canonically corresponds a field v on TM 
taking values in TM (the geometric representative of V), we define: 

Definition 1.5. (Covariant Value). We define the covariant value of the second 
order differential equation D, denoted by , as the field in TM taking values in TM 
corresponding canonically to D — Dq. 

Remark 1.6. The covariant value yields at each point Va € TM the acceleration 
at the point a G M of a particle with a trajectory solution of D and having Va as 
tangent vector at the point a G M: 
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where v is the vector field tangent along the trajectory. 

From this definition, Lemma (11.11) . and the definition of Dq it is straightforward 
to prove: 

Theorem 1.7. The horizontal 1-form corresponding to the second order differential 
equation D is related to its covariant value by 

(10) iD^T2 = -a . 

Definition 1.8. Given a horizontal 1-form a, grada is the field on TM taking values 
in TM such that 

(11) igradaT2 = a . 

In particular, when a = dU for U G C°°(M), grada equals gradU. 
So, Theorem gives 

(12) = -grada . 

Remark 1.9. When T2 is a second order covariant symmetric tensor on M, we can 
define on TM a "Liouville form" ^ by ([1]), and a closed 2-form 002 = dd. Fomulae ([2]), 
(E]), dl]), dSj) remain valid, as also Lemma (11. ip . In Theorem (II. 3p is still true that ([8]) 
assigns a horizontal 1-form a to each second order differential equation Z), although 
not in a one to one way, as a rule. 

1.10. Local coordinate expressions. 

Consider an open set of M with coordinates g^, . . . , and the corresponding open 
set in TM with coordinates g\ . . . , g", g^, . . . , g". 

If the expression in local coordinates of T2 is 

(13) T2 = (7,fc(g) dq^dq^ 

then the local equations for the isomorphism TM ^ T*M are 

(14) Pj = 9jk{q) g'. 
The Liouville form in TM is given by 

(15) 9 = g,k{q) q'dq^ 
and the symplectic form in TM by 



(16) UJ2 = 9ikdq^ A dq'' + ^q^dq^ A dq'' . 

oq'- 



For the kinetic energy we have, locally, 
(17) T = ^g,,q^q\ 
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SO that 

. , dT 
(18) ft = 5^ . 

Let the second order differential equation D be given by 



where the g*'s are given function of g's and g's. The local expression (fT6l) for uj2 gives 

Zd^2 = -gijq^dq' + gijq'dq^ + ^q\q^dq^ - qHq^) . 
On the other hand, we have 

dT = g,,q'dq^ + ^ ^q'q'dq' . 
Replacing this equality into the precedent expression we get 



tnUJ2 = -rfT + 1 ^q'q^dq' + g,,q^dq^ + ^q\q'dq^ - q^dq') 



= -dT + + ^ + ^ - ^ ) ^'^'dq 

= -dT + {gikq' + Tij,kq'q^)dq^ , 

where the F's are the Christoffel symbols of T2. 
Now 

(20) a = -gikiq' + Tyy)dq' 

is the horizontal 1-form related to D by formula (JH]). 
For the geodesic field we have, 



(21) z)o = ,-|--rS,4Y| 

and, finally, the covariant value of D is 

(22) ^^ = (^' + r!,^V)|r. 
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2. Newton-Lagrange mechanics of a free system 

Definition 2.1. (Mechanical System). A mechanical system {M,T2,a) is a man- 
ifold M (the configuration space j provided with a pseudo-Riemannian metric T2 (to 
which corresponds a kinetic energy T by U.^) ). and a 1-form a on TM (the work- 
form or force formj. 

Postulate 2.2. (Newton Law). The tangent field D on TM such that 
(23) iDUJ2 + dT + a = Q 

is a second order differential equation. 

From Theorem (11.31) . we have 

Theorem 2.3. The force-form a of a mechanical system obeying Newton's law is 
horizontal. 

Definition 2.4. (Work). A trajectory of the mechanical system is a curve in TM 
solution of D. The integral of a along a given trajectory is called work done by the 
system. 



Contracting with D in fl23|) . we get 

(24) DT + ioa = 0, 
and, by integrating along a trajectory c, 

(25) j{dT + a) = 

which is an expression for the law of conservation of energy: the work done by the 
system equals the loss of kinetic energy. 

Definition 2.5. (Conservative System). When a is an exact differential form, 
the mechanical system (M, T2,a) is called conservative. The function U (determined 
up to an additive constant) such that a = dU is called the potential energy of the 
system. 

Theorem 2.6. In a conservative system {M,T2,dU), the potential energy U is a 
function on M (lifted to TM). 

Proof. dU is horizontal. □ 

Definition 2.7. (Hamiltonian). In a conservative system {M,T2,dU), the function 
H = T + U in TM is called the total energy or hamiltonian of the system. 



Equation (l23l) for a conservative system is 
(26) iDUJ2 + dH = 
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which, read in T*M, is the set of Hamilton's canonical equations. 

Equation ( |T2i) corresponds to the classical Newton's Law F = m ■ ~a . For a 
conservative system: 

(27) D^ = -gradL/. 

or, in local coordinates, 

.^oN d'^q' , dq' dq^ ,^ dU 

For the geodesic field (a = 0), we obtain the equations of geodesies for T2 

as a generalization of the classical Euler's theorem (see [1]) for the movement of a 
point on a surface in absence of external forces. 

For a general 1-form of force a, the corresponding field D differs of the geodesic 
one, Dg, by the vertical field D — Dg', it is natural to consider this field as the "cause" 
bringing the system to move out the geodesies. So, we give the following 

Definition 2.8. (Force). A force on M is a vertical field V on TM. The field v 
on TM which takes values in TM corresponding canonically to V , will be called the 
geometric expression ofV, or the geometric representative ofV. The force associated 
to a second order differential equation D by the metric T2 is D — Dg- Its geometric 
expression is D"^ . 

3. Constrained systems 

Definition 3.1. (Constraints). A constrained mechanical system is a mechanical 
system (M, T2,a) together with a Pfaff system A on TM. Each 1-form in is a 
constraint . 

Postulate 3.2. (Newton-Lagrange Law). Given the constrained mechanical sys- 
tem {M,T2,a,A), there exists a tangent vector field D on TM, which is a second 
order differential equation, and satisfies the congruence 

(30) i-0^2 + dT + a = mod A 

and also the principle of virtual works (stated below Iji3.4\ )). 

By Theorem (12. 3p a is horizontal. Since D is a second order differential equation, 
the /3's in A transforming (130|) into an equality, have to be horizontal. Thus, we 
can always replace A by its intersection with the space of horizontal forms on TM. 
In the sequel, we will assume that A is a Pfaff system of horizontal forms on TM. 
Therefore, for each j3 in A, it makes sense to consider the function /3 on TM. 
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Definition 3.3. (Admissible State.) An admissible state for the constrained me- 
chanical system (M, T2,a,A) is a point Va € TM such that I3{va) = for each j3 in 
A. The set consisting of all admissible states is denoted by C. 

Principle of virtual works 3.4. A curve in TM which is a solution of D and passes 
through a point Va € C, remains entirely in C 

When £ is a submanifold of TM this is equivalent to D being tangent to C 

Remark 3.5. Congruence fl30l) is a form of D'Alembert principle of equilibrium 
between applied, inertial, and constraint forces. Definition (13.31) selects as admissible 
velocities those for which, in the corresponding "infinitesimal displacements", the 
constraint forces do not work. The principle (13. 4p means that the system remains in 
admissible states: the constraint forces never work. (See, e.g. [0], Sect. 10). 

3.6. Local computation for D. 

We take local coordinates in M and a local basis . . . , (3r} for A in TM: 

(31) A = BkM\ ■■■A^A\-- ■,ql dq' (fc = 1, . . . ,r) 

Each corresponds to a vertical tangent field Vk (a constraint force) by 

(32) iv^uj2 = Pk 
Let D be the field for the free system (M, T2, a): 

^0^2 + dT + a = . 
Then, the congruence ( 130|) can be written as: 

(33) D = D + AVi + ■ ■ ■ + yVr , 

where the A'^'s are certain (local) functions on TM, called Lagrange multipliers. Our 
problem is to compute the A's such that D satisfies the principle of virtual works 

dsai). 

We assume that C is a submanifold of TM. In this situation we compute the A's 
by requiring D to be tangent to C. 

We have the local expressions 

(34) V^, = 4(g^...,g^g^...,g")A (A; = 1, . . . , r). 



The local equations of C are /5i = 0, . . . , /5r = 0. Then, the conditions for D being 
tangent to C are: 

(35) X^a^^!^^ ^ -D{Bi,q^) mod (A, . . . , /3.) 
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Using coordinates, ( l32l) is: 

(36) gjhal = Bkj 
and (!35|) becomes: 

(37) X'af^^0^ ^ -D{g,Hatq') mod (A, ...,/?,). 

The discussion of ( 1371) is carried out according to the supplementary hypothesis in 
each particular case. We start with the most common one: 

Definition 3.7. (Linear Constraints). A is a system of linear constraints when it 
is generated, as a Pfaff system on TM, by a Pfaff system Am on M (lifted to TM). 

Let r be the rank of A^j. For each a G M, the vectors Va G T^M annihilated by 
the forms in {AM)a form a (n — r)-dimensional subspace Ca C TaM. The collection of 
all Cai when a varies in M, is the set C of admisible states. Therefore, in the case of 
linear constraints, C is the vector subbundle of TM corresponding to the distribution 
of vector fields incident with (i.e. annihilated by) Am- 

In this case, the in flMl) can be chosen free from the g's, and the left member in 
f l57|) becomes 



(38) \^a!lgjsal5jh = X^a^g^sa^i = \^{vk,vi) , 

where Vk, vi are the geometric representatives of V^, Vi (see Definition (12.81) ). and ( , ) 
is the scalar product with respect to T2. 

When T2 is positive definite, the matrix with entries {vk,vi) is non degenerate, 
because the f 's are linearly independent. It follows that (!37|) can be solved when we 
replace congruence by equality giving unique solutions for the A's by: 

(39) X'{vk,vi) = -D/3i. 

With these values for the A's, the D in ( |33l) gives: 

(40) DPi = Dl3i + A VfcA = D(3i + \\(3i, v^) = D(3i + \\vu v^) = 

In the open set of TM where we work, the field D satisfies the conditions imposed 
by ( 13. 2p . Also D depends upon the basis . . . ,/9r}, but its restriction to C does 
not depend. 

The local fields D can be used to build, by means of an appropriate partition of 
unity, a global field D on TM satisfying ( 13. 4p . The uniqueness oi D \c follows from 
the local uniqueness. 

We have proved: 
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Theorem 3.8. When T2 is positive definite and A is a linear system of constraints, 
then there exists on TM a second order differential equation D that satisfies the 
Newton- Lagrange law liS. 2\) . The restriction of D as tangent field to the bundle C of 
admisihle states is fixed by the condition settled in Ili3.4\ )- 

Remark 3.9. (About non linear constraints). In the general case of a Pfaff 
system A of horizontal forms on TM, and assuming T2 to be positive definite, the left 
member of fl37|) restricted to the section of TM is (l38l) . Thus, in a neighborhood of 
the section, there exists a unique solution for the A's when in f l371) we replace = by 
=. And, in the analytic case, there exists a global D, with possible singularities out 
of the section (i.e. "for large velocities"). 

The most important case of linear constraints is that of holonomic constraints: 

Definition 3.10. (Holonomic constraints). A system A of constraints is called 
holonomic when it is linear and generated by a completely integrable Pfaff system Am 
on M. 

In this case, let {dBi, . . . ,dBr} be a local basis for Am- The local equations for 
A are Bi = 0, . . . ,Br = 0, and (locally) A is the union of the tangent bundles to 
the submanifolds of M given by equations Bi = hi, . . . , Br = br {bi constants). Let 

C M be one of these submanifolds. In the embedding A^ C M, the metric T2 on 
M specializes as T2n- In the embedding TN C TM, the Liouville form 6 on TM, 
specializes as the Liouville form on TN corresponding to the metric T2N (this is a 
consequence of ([1])). Then, U2 in TM specializes as uj2n in TN, which is also the 
symplectic form corresponding to T2n- 

The field D on TM corresponding to the mechanical system (M, T2, a. A) is tangent 
to TN, because the local equations for TN in TM are: 

Bi = bi, . . . , Br = br, Bi = 0, . . . , Br = 0, 

and, for each Va G TaN we have Dv^{Bk) = Va{Bk) = and Dy^B^ = 0, by 03.41) . 

Now, in the Newton-Lagrange equation 

r 

i-^uj2 + dT + a = ^\'' dBk 

k=l 

everything specializes to TN, and with such a specialization, we obtain 

%^^2Af + dTN + aN = 0. 

Theorem 3.11. When the constraints are holonomic, the field D corresponding to the 
system {M,T2,a,A), tangent to the submanifold C of admisible states, is also tangent 
to the submanifolds TN C C, for each solution N <Z M of the Pfaff system Am- The 
restriction of D to TN is the field corresponding to the free system {N,T2N,CiN)- 
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3.12. Analysis of general linear constraints. 

Let D he a. second order differential equation, j3 a 1-form on M. Working in local 
coordinates we have 

(3 = Bi{q\...,q-)dq' $ = Bi{q\ . . . , q^) q\ 

and 

= B4 + D{B{) q' = (/?, - T^^q'^q''^) + d{Bi)q' 

where is an horizontal 1-form on TM and can be paired by duality with d. 
For the geodesic field, we obtain 

DG$={V^P,d) . 

Let V be the field such that i^T2 = (3. We have 

d) = (V^(^.T2), d) = T2{V^v, d) = \\M d), 
where 11^, is the second fundamental form associated to the field v. 
Putting all together, we obtain 

(41) D/3 = {(3, D^) + Dg/3 = {/3, D^) + ll,{d, d). 

This formula will be used to clarify the behaviour of a mechanical system with 
linear constraints. 

Let (M, T2,a, A) be such a system, let {Pi, . . . , Pr} be a local basis for Am, and 
Vk = grad/^fc, Vfc the vertical representative of Vk, so that ivk^2 = Pk (Lemma (11. ip ). 
We write Ilk for the second fundamental form corresponding to Vk- 

By a suitable choice of the basis, we set 

{vh,Vk) = Shk, {h,k = 1,.. . ,r). 

Then, equations fl39l) for the Lagrange multipliers are 

A' = -DPi, 

when D is the field for the free system (M, T2, a). And the field D for the constrained 
system is, according ( l33l) 

r 

(42) D = D-J2iDPk)Vk, 

k=l 
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From (HTj) and fll2]) . we obtain for the covariant values, 

r r r 

(43) = D^- J2iD$k)vk = D^ - Y.^Vk,^'')vk -Y,^\u{id)vk 

k=l k=l k=l 

The sum — 'Yl^=i{'^ki D^)vk is the orthogonal projection of on the distribu- 
tion C (in this case, the orthogonal distribution to f i, . . . , Vr). The sum — X]fc=i IIfc('^, d)vk 

is the only term surviving when D is the geodesic field (a = 0); then, it is Dq^ the 
field corresponding to the constrained system (M, T2,0,A). 

The formula = — ^^^-^ IIfc((i, (i)ffc, in the particular case of holonomic con- 
straints, express a well known geometric theorem: when iV is a submanifold of M, 
the acceleration in M of a geodesic of N is orthogonal to N . 

Formula (HSj) is a related to classical results; see e.g. Prange 0, pag. 558. 



4. The class of time and the calculus of variations 

As usual, differential 1-forms on M are also considered as 1-forms on TM, by 
means of the natural pull-back. A 1-form a on M defines the function a on TM and, 
in the open set where a does not vanish, it defines the 1-form ^. For each second 
order differential equation D we have (^, D) = 1. So, for two 1-forms a, /3 on M we 
have (| - = in the open set of TM where neither a nor 13 vanish. From the 

definition of the contact system (10. 3p it follows that we have in that open set 

" - ^ o 

— = — modiZ . 

a (3 

Definition 4.1. (Class of Time). Let U be an open set of TM non intersecting 
the zero section. A 1-form r defined on U belong to the class of time when at each 
point Va E U it satisfies 

rva = -j-^ modi2^^ 
for each 1-form a on M such that a{va) 7^ 0. 

The 1-forms belonging to the class of time are horizontal. Therefore we can apply 
them to the field d. The following proposition is obvious 

Proposition 4.2. A 1-form r on an open set U C TM non intersecting the zero 
section belongs to the class of time if and only if it is horizontal and satisfies (r, d) = 1. 
Therefore {t,D)=1 for each second order differential equation D. 

Definition 4.3. (Infinitesimal contact transformation). A tangent field 6 on 
TM is a infinitesimal contact transformation if LgQ C (i.e., for each 1-form a in 
Q, Lga is in Q ). 
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When S generates an uniparametric group of automorphisms of TM, the condition 
Ls^l C Q means that this group transforms solutions of Q into solutions of Q, i.e., 
maps into each other curves which are tangent to second order differential equations 
or vertical fields. 

Theorem 4.4. Let S be an infinitesimal contact transformation on TM, which is 
projectable to M as a tangent field. The following two properties are equivalent: 

(1) 6 leaves invariant the class of time: for each 1-form a on M , in the open set 
ofTM where a ^ 0, we have: 

Ls (^^) = mod Q . 

(2) 6 commutes with d 

d o S = 6 o d 

as derivations from C°°{M) to C°°{TM). 

(Note that d o S does make sense because 6 is projectable to M, thus maps C°°(M) 
into C^{M)). 

Proof. Let D be any second order differential equation and let a be a 1-form on M. 
In the open set of TM where d 7^ 0, we have 

= 5(1) = D) = {Ls (-) , D) + (-, [6, D]) . 

The first term in the last sum is for arbitrary D if and only if 6 satisfies ([T]). The 
second term is for arbitrary a if and only if [6, D] is vertical, and this is equivalent 

to □ 

Definition 4.5. (Infinitesimal Variation). An infinitesimal contact transforma- 
tion that is projectable to M and satisfies the equivalent conditions of Theorem Iji4-4\ l 
is called an infinitesimal variation. 

Theorem 4.6. For each vector field v on M , there exists a unique infinitesimal 
variation 5^ on TM which projects to M as v. 

Proof When 6^ exists, (ED of (g^D gives, for each / e C~(M), 

6 J = d{vf) . 

Thus, 5^ is uniquely determined by v. 

Given a coordinate system g-*^, . . . , g" on an open set U C M, let f = a*(g\ . . . , g")^ 
be the local expression for v. Defining 5^ = + d*(g^, . . . , g")^ we check that 
5^ is an infinitesimal contact transformation by computing the Lie derivatives of the 
forms q^dq^ — q^dq^ (local generators of Vt on TM out of the zero section). Then, from 
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its actual definition, 5^ fulfills of (14. 4p in U . The uniqueness shows that the local 
5yS patch together giving a 5^ globally defined on the whole of TM. □ 

Definition 4.7. (Prolongation of a Field). 5„ is called the prolongation of v to 
TM. 

Remark 4.8. The tangent bundle TM is the space of contact elements Ml in the 
sense of Weil [7] . Formula ([2]) of (14. 4p means that 5^ is the canonical prolongation of 
V from M to M|. Also, 5y is the infinitesimal generator of the prolongation to TM 
of the group (or local group) of automorphisms of M generated by v. 

Remark 4.9. (j2]) of (14.41) is the 'rf o 5 = 5 o c?' in the classical texts of Mechanics; 
e.g. in Sommerfeld [6J, formulae (9) and (9a) of page 175. Our formula ([T]) of (14.41) 
is reminiscent of "keeping fixed the time" in the reasoning from Sommerfeld in pages 
175, 176. Sommerfeld attributes to Euler the formula d5 = 6d. For this reason, we 
shall call ([2]) of (14.41) Euler 's commutation formula. 

The starting point for the applications of variational methods in Mechanics is the 
following theorem (see Prange [5], "Zentralgleichung von Lagrange", in page 531). 

Theorem 4.10. Let M be a manifold provided with a pseudo-Riemannian metric T2, 
9 the corresponding Liouville form on TM, T the kinetic energy. For any second 
order differential equation D and any infinitesimal variation 6, we have 

(44) D{e,6) = 6T-{a,6) , 

where a is the work form corresponding to D by Theorem U.3\) . 

Proof. The Cartan formula for the Lie derivative, with Definition (11.21) for T and 
Theorem (11.31) gives 

Lj^e = iode + d{e, D) = i^t^a + 2dT = -dT - a + 2dT = dT - a , 

thus, 

D{9, 6) = {dT -a, 5) + {9, [D, S]) = ST - {a, S) , 
because [D, 6] is vertical as follows from Euler's commutation formula. □ 

Definition 4.11. (Lagrangian function). In a conservative system {M,T2,dU), 
the function L = T — U is called the Lagrangian function of the system. 

Applying Theorem (14.101) to a conservative system, we obtain: 

Theorem 4.12. (Hamilton's Principle). Let {M,T2,dU) be a conservative me- 
chanical system, D the corresponding second order differential equation and L the 
lagrangian function. For each infinitesimal variation 6, we have 



(45) 



D{9,6) = 6L 
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The classical integral version of Hamilton's principle follows from (14.121) as we show 
now. 

Let c: [to, ti] ^ M be a parameterized curve such that its canonical lift to TM is a 
solution of D. Let f be a tangent field on M vanishing at the points c{to) and c{ti); 
let be the prolongation of v. 

Along the lifting of c to TM, D is Then, by integrating fHS]) and using that 
{9,Sy) = {9,v) vanishes at to and ti, we find 

(46) [ {6yL) dt = , 

Jto 

which is the classical integral form of Hamilton's principle. 

Newton equation fl27j) is the "Euler-Lagrange" system for the variational principle 

(USD. 

Remarks 4.13. (About constrained systems). Let (M, T2,a,A) be a me- 
chanical system with linear constraints (13. 7p . Let C C TM be the linear bundle of 
admissible states, (13. 3p . A vector field t> on M is an admissible virtual displacement 
when Va & C for each a G M; then, the infinitesimal variation 6y is an admissible 
infinitesimal variation. 

Such a 6v is not tangent to £ in general. The local uniparametric group generated 
by 6y may transform curves in £ into curves out of C, i.e. kinematically possible paths 
into kinematically impossible paths. 

Let D be the second order differential equation corresponding to (M, T2,a,A). 
According (14.21) . there exists a 1-form /5 in A such that 

i-QUJ2 + dT + a + P = . 

Then (1441) gives 

D{e,5) =5T-{a,S) - 

for each infinitesimal variation 6. When 6 = 6y for an admissible virtual displacement 
V, we have 6) = {(3, v) = 0. When, on top of the above, the system is also 
conservative {a = dU), we obtain 

(47) D{e, 5) = 5,L , 
and by the same argument given for (H6|) . we have 

(48) [\6yL)dt = 

when we integrate along a curve c solution of D and v is any admissible infinitesimal 
displacement null at c(to), c(ti). 

But now, we cannot in general reach all the kinematically possible paths close to 
the given c by means of deformations along solutions of admissible (5„'s. So, we cannot 
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derive the extremality of real trajectories among all kinematically possible paths. See 
e.g. Whittaker [S], pag. 250, for a discussion. 

In the case of holonomic constraints, A admits local basis of exact 1-forms {dBi, 
. . . ,dBr}. When v is an admissible virtual displacement, we have v{Bk) = and, 
from Euler commutation formula, S^Bk = 0. Then, the local equations for C in 
TM: Bk = const, Bk = {k = 1, . . . ,r) are preserved by 6jj- So each admissible 
infinitesimal variation is tangent to C; and, indeed, tangent to the tangent bundle 
TN for each one of the solutions of the distribution C. In this way, Hamilton's 
variational principle specializes, like the Lagrange equations, in the case of holonomic 
constraints. 



5. Time constraints 

Definition 5.1. (Time constraint). A time constraint for M is a connected closed 
hypersurface lA of TM, non intersecting the 0-section, projecting regularly onto M 
and such that the class of time is locally represented in U by exact differential forms. 

Remarks 5.2. The last condition means that for each f ^ G W there exists an open 
neighborhood Ui of Va in U and an exact 1-form ri on Ui that belongs to the special- 
ization to Ui of the class of time. 

Let / G C°°{Ui) be such that ri = df. Since the forms in the class of time are 
horizontal, so is df. Therefore, by restricting Ui to a smaller neighborhood of Va if 
necessary, / is the lift to Ui of a function on M. So the local representatives of the 
class of time in U are differentials of functions on M (lifted to TM). 

Lemma 5.3. Let U be a locally closed submanifold of TM, dimU = n + r (r > 
1), projecting regularly onto an open set of M. For each Va & U there exists a 
neighborhood of Va in TM and r linearly independent tangent vector fields on this 
neighborhood, which are second order differential equations tangent to U. 

Proof. Let (g^, . . . , g", g^, . . . , g") be local coordinates on a neighborhood of Va in TM. 
The projection of W to M being regular, and reordering the coordinates if necessary, 
we can find local equations for U of the form: 

cf^' = F^^\q\...,q\q\...Al (A; = 1, . . . , n - r) . 
Then, the fields 

satisfie the requested conditions. □ 
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Lemma 5.4. Let U he as in ( 15.31) . and suppose that U does not intersect ^-section. 
Given f G C°°{M), the differential df represents in lA the class of time if and only if 
/ = 1 identically on U. 

Proof. When / = 1 identically on U, we have y = df , so df represents the class of 
time. 

Conversely, let df represent the class of time in U. From (15.31) . for Va & U there 
exists a second order differential equation D tangent to W in a neighborhood of Va 
in TM. Because df belongs to the class of time in U, we have {df, = 1 in a 
neighborhood of Va in U. In particular, we get {df,Va) = 1 and /(fa) = 1- D 

Lemma 5.5. Assume n > 2. Let U be a time constraint for M. Then, for Va E U 
there exists a neighborhood U{va) in TM and f G C°°{M) such that the equation for 
UnU{va) tn Uiva) ts f = l. 

Proof. From (15.21) . there exists / G C°°(M) such that df represents the class of time in 
a neighborhood of Va in U. From (15. 4p it follows that / = 1 in this neighborhood. We 
have d^^f ^ because daf 7^ 0. Thus / = 1 is the local equation for an hypersurface 
of TM in a neighborhood of Va- This hypersurface has to coincide with U in a 
neighborhood of f^. □ 

Lemma 5.6. With the same conditions as in Ii5. 5\) . let g G C°°{M) be such that 
g = 1 on a neighborhood of Va in U. Then f and g differ by an additive constant on 
a neighborhood of a in M. 

Proof. Let g^, . . . , be local coordinates for M on a neighborhood of a, taken such 
that f = q^- Then, in the neighborhood of Va in TM, where the local equation of U 
is q^ = 1,U consists of the tangent vectors to M of the form 

The equation g = 1 onU {va) gives 

for arbitrary A-' G M and 6 in a neighborhood of a. This implies g = q^ + const in this 
neighborhood. □ 

Lemma 5.7. Let U be a time constraint for M , and n > 2. For each a G M, the 
fiber Ua = U n TaM is an affine hypersurface ofTaM. 
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Proof. As W is a closed hypersurface of TM that projects regularly onto M, each fiber 
Ua is a closed hypersurface of T^M. From Lemma flS.Sp it follows that Ua is locally 



Theorem 5.8. There exists a canonical one to one correspondence between closed 
1-forms without zeros on M and time constraints for M. The 1-form r defines the 
time constraint U with equation f = 1 in TM, and for each time constraint U there 
is a unique r which defines lA by r = 1. 

Proof. First suppose that n = dimM > 2. Let U he a time constraint for M. Lemmas 
f l5.5p and (15 .7^ show that for each a G M there exists a neighborhood Mi of a in M 
and a function / G C°°{M) such that the equation of W fl TMi in TMi is / = 1. 
Lemma (15.61) shows that df is uniquely determined by this condition (reducing Mi if 
necessary). Patching together the d/'s we obtain a closed 1-form r on M such that 
U is defined in TM by the equation f = 1. And this r is unique. Conversely, given 
T closed, without zeros on M, the equation f = 1 defines a time constraint for M. 

Now, let dim M be 1. Then, a closed hypersurface U of TM which projects regularly 
onto M and does not intersect the 0-section is the same that a vector field v without 
zeros on M. The 1-form r defined by (r, v) = 1 defines U hj t = 1. □ 

Remarks 5.9. According to Theorem (15.81) . when M is compact and simply con- 
nected, there are no time constraints for M, because a closed 1-form on M is exact, 
and necessarily has some zeros. In that case, time constraints can be considered for 
open submanifolds of M and non-vanishing exact forms therein. 

When M is compact, there is no "time function", i.e., a function t G C°°(M) such 
that dt defines a time constraint, because dt has always zeros. So, if M is compact, 
"time has always periods". 

6. Modification caused by a time constraint in a free system 

Let (M, T2, a) be a free mechanical system, and D the corresponding second order 
differential equation, so that 



Let U C TM be the time constraint defined by f = 1, where r is a closed 1-form 
without zeros on M. 

Let gradr be the tangent field on M such that igradr72 = t (11.81) . and Gradr the 
vertical field canonically associated to gradr, so that (Lemma (11.11) ) iGradT^^2 = t- 

Theorem 6.1. When gradr is non isotropic for T2 at any point of M (in particu- 
lar when T2 is positive definite), then there exists a tangent vector field D on TM, 
satisfying the conditions 



affine, thus affine. 



□ 



(49) 



iD^2 + dT + a = 



THE STRUCTURE OF TIME AND INERTIAL FORCES IN LAGRANGIAN MECHANICS 23 



(1) i-j^2 + dT + a = Q modij) 

(2) D is tangent to the hypersurface lA. 

Such field is a second order differential equation and its restriction as tangent field to 
U is uniquely defined by the conditions (Qp, 

Proof. Let D be defined by 

(50) D = D-- -Gradr, 

(grad r, grad r) 

where (grad r, grad r) is the scalar product defined by T2. Inmediately D satisfies ([1]), 
because D satisfies f H9l) . 

Also D satisfies ([2]) because we have 

(Grad r)(f) = (grad r,r) = (grad r, grad r) , 

thus = 0. 

Like D is a second order differential equation, because Gradr is vertical. 

Any tangent field on TM that satisfies ([1]) should have the form D + AGradr 
(A G C°°(TM)), and the condition ([2]) determines uniquely the value of A on U. That 
proves the uniqueness of D on W. □ 

Indeed the same proof shows: 

Theorem 6.2. When gradr is non isotropic for T2 at any point of M (in particular 
when T2 is positive definite), there exists a tangent vector field D on TM, which 
satisfies the conditions 

(1) i-QUJ2 + dT + a = raod^r) 

(2) D is tangent to each hypersurface f = c (c eM.) ofTM. 

Such a field is a second order differential equation and is given by formula (E^. 

Remark 6.3. According to (16.2p . the vector field D, when restricted to the hyper- 
surface f = c {c E M.) gives the evolution equations of the mechanical system with 
time constraint defined by ^ r, when c 7^ 0, or holonomic constraint f = in the limit 
case c = 0. In this sense, ordinary holonomic constraints result from "freezing the 
evolution of the system with respect to r" . 

An elementary example. Let M = W^, T2 = dx^ + dy'^ + dz'^, a = and the time 
constraint r = 1 (r = a/o;^ + + z'^). 

Physical Interpretation: a moving point is constrained to be at a distance r = t from 
the origin, at each instant t, with no other external force acting upon it. 
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We have 



X d yd z d 
r dx r dy r dz 



gradr = —— — I — — — I — — , ||gradr|| = 1 
Thus, the geodesic field, modified by the time constraint is 
D = D — (Df) Gradr, where D 



d d 
dy 



dx 



.d_ 
dz 



(the geodesic field for the Euclidean metric on 
A direct calculation gives 



D = D 



[xx + yy + zz) 



x^ + y'^ + z^ 



d 



d 



d 



dx dy dz 



D is tangent to the manifolds r = const in TR^. 

i) On f = (ordinary holonomic constraint r = const), we have 



Dr 



. d . d . d x'^ + y^ + z"^ 

dx ^ dy dz 



d d d 
dx dy 



dz ^ 

sphere of radius 



which is the geodesic field when restricted to each TS C TM^, S 
tq with the center at the origin. 

In the submanifold r = of TM^, to which Dq is tangent, the functions r and 
V = a/A^ + + are first integrals for Dq, and the system of differential equations 
defined by Dq is 

(Px 
1^ 



- I - I a; 



(fy (Pz (V 



which correspond to the classical centripetal force which generate the motion along 
geodesies of the sphere. 

ii) On r = 1 we have 



- d . d . d l-(x^ + f + z^) 

D,=x— + y— + z— + ^ ^ ^ 

dx dy dz 



d 



d 



d 



dx dy dz 



that corresponds to a force depending on position and velocity, which points to or 
out of the origin according to the value of the velocity. 

Remark 6.4. Later on, we show that, in local coordinates, fISUl) and f = 1 give the 
differential equations for a mechanical system with kinetic energy and force depending 
on time, as they are found in the classical literature, e.g. Prange [5], page 556. These 
equations are usually derived from the "Zentralgleichung" in (jH]) when applied to 
arbitrary 5's that keep fixed the time, i.e. such that (r, 6) = 0. But, apparently, there 
appears to be no mathematical justification for this choice of 5's, although some 
authors try to give a physical one; e.g. Sommerfeld [6], page 65, Nordheim [1], pages 
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48, 49. In that follows, we give a precise mathematical justification, starting from the 
considerartion of time as a constraint. 

Lemma 6.5. The prolongation 6y of a vector field v on M to TM is tangent to a 
time constraint lA C TM defined by f = 1 if and only if (r, v) is a constant function 
on M. 

Proof. The problem is local in M, so we can take t = df for some / G C°°{M). 
Euler commutation formula gives 

(51) 6yf = 6yf = d{vf) = d{T, v) . 
So, 5v tangent to W <(=^ (i(r, t>) = on lA. 

Thus when (r, v) is a constant function on M, 5^ annihilates f, and is tangent to 

U. 

Conversely, when 5^ is tangent to U, by taking in fl5T]) the value of d at each Wa G W, 
we obtain Wa{T,v) = 0. But, for each a G M, the set consisting oi Wa & U is the 
affine submanifold of T^M defined by {wa^, t) = 1, and such submanifold generates 
the whole of T^M. Therefore, (r, v) is constant on M. □ 

Definition 6.6. (Admissible infinitesimal displacement). A vector field v on 
M is an admissible infinitesimal displacement for the time constraint U when 5^ is 
tangent to U. In that case, 6^ is called an admissible infinitesimal variation. 

From Lemma (16.51) it follows 

Theorem 6.7. Let (v) be the C°°{M)-module generated by the vector field v on M. 
The necessary and sufficient condition for every field in (v) to be an admissible infin- 
itesimal displacement for the time constraint U (f = 1) is {t,v) = 0. 

For variational problems with fixed end points, the admissible infinitesimal dis- 
placements has to admit multiplication by arbitrary functions on M. Thus, Theorem 
(16. 7p offers a precise mathematical justification for the procedure of "keeping fixed 
the time" . Time must be considered as a constraint. 

6.8. Local expressions 

Let us take local coordinates on M, (g°, q^, . . . , g™) (n = m+1) with time constraint 
g° = 1; so, g° is time. We have gradr = gradg°, and Hgradrp = g^'^ (coefficient of 
the dual metric). The differential equations for (M, T2,a) with no time constraint, 
with a = Aidq^, are 

(52) g,^q^ + F^^^.g^^g^ + A = (z = 0, 1, . . . , m) . 
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When we impose the constraint = 1, the field in (l52l) changes to D, given by 
(150|1 . to which corresponds a force form a that satisfies 

a = a mod (dq^) . 

Therefore changing D hj D does not modify Equations / f^) for indexes i ^ (the 
"spatial indexes"). By separating spatial indexes /x = 1, . . . , m from the temporal one 
0, and putting g° = 1, we obtain the differential equations for the system with g° as 
time: 



(53) 



9^,^q'' + T.^^^q'^ct + 2T^Q^^cf' + Too,^ + = , 
see e.g. Prange page 556, Eq. (64. a). 

The "contravariant " form of the equations results from writing (150|) in coordinates. 
The equations are 

'cP = 1 

q^ + T^^^q'^cf + 2T%qP + T^q + (^^^ + 12! ^o^) = g . 



(54) 



Notice that the difference with the equations given by Prange [5], page 556, Eq. 
(64. b) is only formal. In Prange's equations the translation from subindices to su- 
perindices is done by using only the "spatial" matrix {gfj,u) instead of the whole matrix 
{gij) that we use. 

When we use as time instead of g°, the equations in the constrained manifold 



g° = c are 



(55) 



q^ = c 



g^u'(r + ra.,^g"r + 2cr,o,Mg" + Poo,/. + A^ = Q. 



When c = 1, these give (153|) . And when c = 0, (155!) are the equations for the 
holonomic constrained system (M, T2,a,g° = const). 

7. Linear constraints depending on time 

Definition 7.1. (Linear Constraints Depending on Time). Let {M,T2,a) be 
a mechanical system. A system of linear constraints depending on time is a Pfaff 
system Am on M together with a closed 1-form t, without zeros on M, and such that, 
for each a e M , Ta ^ {^M)a- 

A mechanical system with this type of constraints is denoted by (M, T2, a. Am, t). 

Definition 7.2. (Admissible State). An admissible state for the constrained me- 
chanical system (M, T2, a, Am, t) is a point Va G TM such that, for each Pa € (AM)a, 
we have {Pa,Va) = 0, and {Ta,Va) = 1. 
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The set consisting of all the admissible states is a submanifold V C TM, V = £nW, 
where L is the vector distribution annihilated by Km-, and U is the time constraint 
defined by f = 1. For each a G M, the fiber Va is the affine hypersurface of defined 
by equation (r^, w^) = 1. 

Theorem 7.3. (Newton-Lagrange Law). Let (M, T2, a, Am, t) he a mechanical 
system with constraints depending on time, and let T2 be positive definite. There exists 
a vector field D on TM, which is a second order differential equation, satisfying the 
congruence 

(56) iT)0J2 + dT + a = Q mod {Km, t) 
and the Principle of Virtual Works; 

(57) D is tangent to the manifold V of admissible states. 

The restriction of D to V is uniquely defined by these conditions. 

Proof. Denote by D the second order differential equation corresponding to the free 
system (M, T2,a) such that 

(58) iDUJ2 + dT + a = . 

Let U he a coordinate open set in M such that Am admits on f/ a basis {Pi, . . . , Pr}', 
denote r = /3o, and Vk = grad/3fc. Let Vk be the vertical representative of Vk, so that 
ivk^2 = Pk (A; = 0, 1, . . . , r). Similarly as in Theorem ( 13.8p . we can find functions 
on TM (the "Lagrange multipliers" ) such that the vector field D defined on TU by 

r 

(59) D = D + ^ AVfc . 

A:=0 

satisfies Dpo = Dpi = ■■■ = DPr = on f/. 

D is a second order differential equation and satisfies both (1561) and (1571) on U. 
The field D built in this way depends on the choice of the local basis {Pi, . . . , Pr} for 
Am, but its restriction to the submanifold V HTU does not. By using an appropriate 
partition of the unity in M (lifted to TM) we find a D fulfilling fl56|) and fl571) and 
uniquely defined on V by these conditions. □ 



7.4. Holonomic constraints depending on time. 

Suppose Am to be completely integrable. Let U be as in the proof of (17. 3p and 
Pk = dBk {k = 0,l,...,r). 

The equations for V fl TU in TU are 
(60) Bo = l,Bi = ■■■ = Br = 
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The field D in (17.31) satisfies 

(61) IDBq = DBi = ■■■ = DBr = onVnTU. 

Let N C U a. solution of A^, given by equations 

(62) Bi = bi,...,Br = br ibkER). 

The equations for V fl TN are ( l62i) and ( l60l) . Applying D to equations fl^ we 
obtain equalities on V fl TN because of ( 1601) : and applying D to ( !60l) also yields 
equalities on V fl TA^ because of fl6Tl) . Therefore, D is tangent to V fl TN. 

For this reason, in the Newton-Lagrange equation fl56l) we can specialize to the 
submanifold V fl TA^ and obtain 

(63) hN^2N + dTN + aN = mod (tn) 
and 

(64) Dttv = . 

Equations (1^ and show that is the second order differential equation for 
the mechanical system {N,T2N,C(]\f) with time constraint = 1- We have proved 
the following: 

Theorem 7.5. Let {M,T2,a) be a mechanical system with T2 positive definite. Let 
Am be a completely integrable Pfaff system on M. Let r be a closed 1-form without 
zeros on M, f = 1 the time constraint defined by r. Let D be the second order 
differential equation corresponding to the system (M, T2, a, A^/, r). Then, for each 
submanifold N G M solution of Am, D is tangent to the time constraint f^v = 1 
in TN, and the restriction of D to this submanifold is the second order differential 
equation corresponding to {N,T2N,CiN,'^N) by the Newton-Lagrange law Ii6. 

Remark 7.6. Theorem (17.51) means that time constraints specialize properly in the 
presence of holonomic constraints: we can introduce the time constraint after spe- 
cializing T2 and a to the submanifolds solutions of A^. 

Remark 7.7. It is clear that the precedent methods can be applied in some cases 
in which T2 is not positive definite, under appropriate hypothesis (e.g. when gradr 
is a field of the distribution C and the restrictions of T2 to C and to the orthogonal 
complement of C are positive or negative definite). 



8. Reference frames and inertial forces 



8.1. General principles. 
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In the context of Lagrangian Mechanics it is natural to consider as reference frame 
of a manifold Ai a manifold isomorphism ip: TZ —>■ Ai that transports mechanical 
structures from A4 to TZ, where we "read them". When TZ and Ai are provided 
with pseudo-Riemannian metrics, the geodesic field of A4, once transported to TZ, 
is a second order differential equation that, by substraction of the geodesic field of 
TZ, gives a vertical field on TTZ: the inertial force caused by ip. The most relevant 
dynamical elements appear when TZ and Ai are endowed with time constraints which 
correspond to each other by if. 

Uniparametric automorphisms groups of a given manifold M are a particular case 
when we consider one such group as an automorphism of the manifold = M x 
M. Classical examples are uniparametric groups of isometrics of M'^, which generate 
centrifugal and Coriolis forces. The same method allow us to deal with more general 
groups, like the one of dilatations that are apparently ignored in the literature. 

Example 8.2. (Inertial forces caused by automorphism groups) 

Let if : TZ ^ Ai he a,Ti isomorphism of manifolds. Let T2 be a metric on Ai, and 9 
the correspondig Liouville form (see (11.91) ). It is easy to prove that ip*9 is the Liouville 
form for ip*T2 on TZ; thus, if U2 = dO, ip*uj2 is the 2-form associated to ip*T2 on TZ. 

In particular, let T2 be non degenerate (i.e., a pseudo-Riemannian metric) and ip 
be an isomorphism. Then, ip^, applies the second order differential equation D into 
D, related by: 

ini^* UJ2) + d{^*T) + ^* {a) = 

when 

%i^2 + dT + a = . 

When a = 0, we observe that the geodesic field for ip*T2 on TZ is transformed into 
the geodesic field for T2 on A^. When TZ is also equipped with a metric, the vertical 
tangent field on TTZ, difference between the geodesic field for ip*T2 and the one of the 
metric given on TZ, is the inertial force produced by ip. 

Let us consider the case 7^ = = M x M, where M is an n-dimensional manifold. 
And let M X M — i> M X M an uniparametric automorphism group: 

ip{t,a) = {t,ipt{a)) , 

where, for each t eM., ipt: M ^ M and ipo = Id, ipfO ip^ = ipt+si as usual. Denote by 
u the infinitesimal generator of ip. 

Given a metric T2 on M, we endow M x M with the metric T2 = dt^ © T2. We give 
a formula for ip*T2. 
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We have </'*(^) = ^ + (we use the product structure on M x M to extend u 
from M to M X M, and f from M to M x M). Thus 



•fit (a) ■ 



For a field v on M, 



9 

(^«a(¥'t^2),^^a) = {dt ® iua{^*tT2)) ( 737, ^^a 



(t,ipt(a)) 

dt 

For two fields f , on M: 

\ / {t,a) \ /a 

We obtain: 

(65) = (1 + rft' + 2rft^„„ (y,:T2) + (^:T2), 



To illustrate with simple numerical examples, let us take M = R , T2 = dx + dy 
then T2 = dt"^ + dx^ + dy"^. And the most classical groups: 

Translations: ip\ {t,x,y) ^ {t,x + t,y). 
In this case, we have u = d/dx and 

(66) ip*f2 = 2dt^ + 2dtdx + dx^ + dy'^ . 
Rotations: ip: {t,x,y) {t,x cos t — y sin t, x sin t + y cos t) . 

^1 ^ II l|2 2,2 *Tn rri 

'^^ ~^dx dy ' +1/ , <ftT2 = T2, 

thus, from ( l65l) . we obtain 

(67) (p*f2 = {l + x^ + y^) dt^ + 2 dt{-ydx + xdy) + dx^ + dy"^ , 
that can also be calculated directly. 

Dilatations: ip: {t,x,y) i— > (t, e*x, e*?/), 

« = a;^ + 2/^, \\uf = x^ + y\ ip;T2 = e^'T2, 
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l|w||Jt(a;,2/) = e^*(^'^ + y^) > ^u{VtT2) = e^\xdx + ydy) , 

thus 

(68) ip*f2 = [1 + e^\x^ + y'^)] dt^ + 2 e^* + ydy) + e2*((ix2 + dy"^) . 

We can use these formulas in order to calculate the inertial forces associated to 

each of these groups. In each case, we consider in M x M the time constraint t = 1 
(M is the clock for the system). The geodesic field for the given metric T2 is Vq ~ 
td/dt + xd/dx + yd/dy, which is tangent to the given constraint. 

Translations. The transformed metric (p*T2 has constant coefficients. Then all 
Christoffel symbols are zero, and the geodesic field for ^p*T2 is the same as the one 
for T2. The inertial force is 0. 

Rotations. A routine tedious computation for the field V such that 

iv{^*^2) + d{v*f) = Q 
{V transforms by ip^ to the geodesic field on T(M x M^)) gives 

^ = + -^TT + + + 2 iy)^ + (y - 2 tx)-^ 

ot ox oy ox oy 

which, in the locus of the time constraint t = 1 is 

The difference with the geodesic field is the inertial force corresponding to the group 
of rotations. The differential equations for this force are 

d'^x dj^v 

(70) ^ = - + 2*. J = »-2*. 

where we observe the sum of a centrifugal and a Coriolis force. 

Dilatations. A similar computation gives for the field T) which transforms by 
into the geodesic field, when restricted to time constraint i = 1: 

d d 

(71) V\i^,^VG-{x + 2x)—-iy + 2y)-. 

To which corresponds the differential equations for the inertial force caused by </?: 

d'^x d'^v 

(72) _ = _(, + 2x). ^=-(v^2y). 

We do not give a physical interpretation of this formal result. 
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Definition 8.3. (Inertial frames). Let TZ, M. he pseudo-Riemannian manifolds. 
A reference frame ip: TZ ^ M. is called inertial when the force caused by ip is zero. 
This applies in particular to uniparametric automorphism groups. 

Equivalent characterizations of inertial frames are: (p is inertial if and only if (p^, 
transforms the geodesic field of TTZ into the geodesic field of TAi. And also, if and 
only if it transforms one to another the Levi-Civita connections on TZ and Ai. 

Example 8.4. (Inertial groups on M"). 

The metric we consider in M" is the usual Euclidean one. We prove that the only 
inertial uniparametric groups on are the translation ones. 

Retaining the notations of (18.21) . for (p to be inertial, T2 and (p*T2 must define the 
same Levi-Civita connection on M x M". Thus, when written in cartesian coordinates, 
all the Christoffel symbols for (p*T2 must to be zero; and then, the gij for ip*T2 are 
constants. Looking at the last term in (l65l) . this gives <plT2 = ^qT2 = T2. Then, since 
the coefiicient in the second term in fl65|) are constant, we obtain that u has constant 
coefiicients. We conclude that ip is the translation group generated in by u, q.e.d. 

8.5. Preservation of the equations of motion. 

In the physical literature, inertial frames are frequently referred as those which 
"preserve the form of the equations of motion", without defining a precise meaning 
to "preserve" . We treat this question in this section. 

Let (7^, T2) and {A4,T2) be pseudo-Riemannian manifolds, and (p: TZ ^ A4 a 
reference frame (= manifold isomorphism). Let D be a second order differential 
equation on TAi, and a its canonically associated 1-form: 



There are two different natural ways to define a second order differential equation 
on TTZ corresponding hj ip to D: 

a) The second order differential equation D on TTZ canonically associated to 'p*{a): 



(73) 



i-jyuj2 + dT + a = 



(74) 



iDUJ2 + dT + ip*{a) = . 



b) The field Di on TTZ such that (p^Di = D: 



(75) 
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Classically, Di is the field obtained from D by change of coordinates. 

The difference Di — D is a vertical field on TTZ, i.e., a force; when we think of Di as 
D "read in the coordinates of TZ" , Di — D is the force caused by (f in the mechanical 
system {A4,T2,a). The frame may be inertial in the sense of Definition (18. 3p . but 
causes a non zero force when a 7^ 0, as shown in the following 

Proposition 8.6. Let (7?., T2) and {Ai,T2) be pseudo-Riemannian manifolds, and 
tp: TZ ^ M. a reference frame. The force caused by ip in (A^,T2,a) is zero for all 
choices of a, if and only if cp is an isometry. 

Proof. The "if" part is obvious. For the converse, let us take a = 0; then, for the 
geodesic field Dq on TTZ, we have 

(76) iDaUJ2 + dT = 0, iDa{v*^2) + d{p*T) = . 

Then, let a be any horizontal 1-form on TAi, and a = 'p*{a)] let D be given by 
fITil) . and define the vertical field V hj D = Dg + V,so that, from flTillTHl) . 

ivUJ2 + a = . 

By taking the field Di in ( |75l) . the hypothesis Di = D gives, from (17511761) . 

iv{p*uj2) + a = . 

Thus ivuj2 = iv{'^*^2)] then, the geometric representative v oiV satisfies 

ivT2 = iv{p*T2) . 

Since V is arbitrary, this implies T2 = p*T2. □ 

Thus it is not true that the uniparametric groups of translations considered in 
(18. 2p do not produce forces; they do, but these forces may be compensated by time 
constraints, as we will show. 

Definition 8.7. (Preservation of the equations of motion). Let (7^, T2) and 
{A4,T2) be pseudo-Riemannian manifolds with respective time constraints f = 1 and 
r = 1 . Let ip: TZ ^ M. be a reference frame. We say that p preserves the equations 
of motion when (P*{t) = r and, for each second order differential equation D on TAi 
such that Df = 0, and for the corresponding a by (7^ ), and for the field Di on TTZ 
such that pifDi = D , we have 

(77) iDi^2 + dT + V5*(a) = mod^r) . 

Di is tangent to the time constraint f = 1, and the congruence (!77|l means that 
Di is the modification of D in (!74|) caused by the time constraint, in the sense of 
Theorem (16.11) . 
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Theorem 8.8. An uniparametric automorphism group ip: RxM^MxM preserves 
the equations of motion if and only if it is inertial and also each ipt: M ^ M is an 
isometry. 

Proof. The problem reduces to prove that for an inertial if, the condition of preserving 
the equations of motion is equivalent to be an isometry group of M. 

For any uniparametric automorphism group ip the geodesic field Dq satisfies Dot = 

0. 

Let D = Dg + V he a second order differential equation on T(R x M) such that 
Dt = 0, i.e. Vi = 0. We keep notations from (18. 2p and define 5 by 

ij)UJ2 + dT + a = 0, or iyUJ2 + 5 = 0. 

Let Di be the field defined by p^^Di = D, and put Di = Dg + Vi. As is inertial, 
we have p^Vi = V, thus 

(78) iv,{p*U2)+p*ia) = . 

Let D = Dg + V he defined by 

(79) iz)^2 + c?T + </?*(5) = , 
thus 

(80) iv^2 + P*ia) =0 
and from (178|) . 

(81) ivu}2 = ivi{'f*(^2) ■ 

The condition for (p to preserve the equations of motion is 

(82) iD,{^2) + df + ip*{a) = mod (dt) 

Substracting (|79l) from (|82l) . this condition is 

ivii^2) - = mod ((it) 

or, from flHTl) 

ivi{i02) = Vi(v5*t^2) mod ((it) . 
By taking geometric representatives, this condition is 

(83) v(T2) = v(v9*T2) mod ((it) . 

Thus the necessary and sufficient condition for (f to preserve the equations of motion 
is that ( l83l) holds for each tangent field vi on T(]R x M) taking values in T(R x M) 
and satisfying vit = (this is Vit = 0). 

Looking at (165!) we see that this condition for ip is equivalent to have PIT2 = T2 
for each t G M. □ 
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Corollary 8.9. The uniparametric automorphism groups of M" which preserve the 
equations of motion are only the inertial groups, i.e., the translation uniparametric 
groups. 
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